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The unoccupied states of complex materials are difficult to measure, yet play a key role in deter- 
mining their emergent properties. We propose a novel technique that can measure the unoccupied 
states, time-resolved Compton scattering, which measures the time-dependent momentum distri- 
bution (TDMD). Using a non-equilibrium Keldysh formalism, we study the TDMD for electrons 
coupled to a lattice in a pump-probe setup. We find a direct relation between temporal oscillations 
in the TDMD and the underlying unoccupied states, suggesting that both can be measured by 
time-resolved Compton scattering. 



PACS numbers: 78.47.je, 79.60.-i, 63.20.kd 

In understanding the emergent properties of complex 
materials, it is insufficient to limit ones' study solely to 
the occupied electronic states. The unoccupied states 
play an important role in determining, for example, the 
nature of the gap in charge- and spin-density waves, the 
absorption properties of semiconductors, or the magnetic 
properties of Mott systems. For example, without knowl- 
edge of the unoccupied states, it is difficult to measure 
both the full width as well as the q of a gap in the elec- 
tronic spectrum. Unfortunately, the most direct exper- 
imental measurements of momentum-resolved states in 
quantum materials, angle-resolved photoemission spec- 
troscopy (ARPES), only determines the occupied states. 
This leaves the alternative methods of inverse photoemis- 
sion spectroscopy (IPES), which lacks the signal strength 
to measure with sufficient resolution, or indirectly infer- 
ring the unoccupied dispersion from two-particle quan- 
tities such as optical spectroscopy. Pump-probe spec- 
troscopy, and in particular time-resolved (tr-) ARPES, 
can measure part of the unoccupied states, but is limited 
by the inherent competition between temporal and en- 
ergy resolution, and cannot map out the entire Brillouin 
zone in most cases. 

Here we consider a different quantity, where one sacri- 
fices the energy information in favor of high time resolu- 
tion, while maintaining high momentum resolution. We 
propose to measure the (gauge-invariant) time-dependent 
momentum distribution (TDMD) of quasiparticles, de- 
fined by (cj c (t)ck(t)) = rik(t), which is similar to the so- 
called Wigner distribution. It contains a rich amount of 
information about the behavior of pumped quasiparticles 
while avoiding some of the resolution-related complica- 
tions of other measurements. 

The TDMD is commonly measured via time-of-flight 



absorption images in cold atomic gases, where it was re- 
cently used to map out the Fermi surface of a gas of 40 K 
atoms. pQ In the context of condensed matter systems, 
we propose two new types of measurements to access 
the TDMD. First, we propose to extend the technique 
of Compton scattering for solids into the time domain [2]. 
Compton scattering has a long history as a probe of 
the electron momentum distribution for solids in equilib- 
rium [3]-[5]. With the advent of x-ray free-electron laser 
sources, the photon energy of the ultrashort x-ray probe 
pulses has been extended into the hard x-ray regime, 
where the fractional Compton cross-section becomes ap- 
preciable. The planned upgrade of the Linac Coherent 
Light Source at SLAC National Accelerator Laboratory 
will extend the energy of the x-ray free electron laser to 
beyond 20 keV. There is also significant progress in crys- 
tal optics at these energies, such that time-resolved high 
resolution Compton scattering with momentum resolu- 
tion of a few percent of the Brillouin zone will be feasible. 
Additionally, recent developments in photonics suggest 
the future availability of ultrashort pulses ranging from 
x-rays to gamma rays [6]. Time-resolved Compton scat- 
tering has the advantage of using high-momentum pho- 
tons, which can readily access the Brillouin zone edges in 
all directions. Secondly, one can access the TDMD with 
tr- ARPES by using similar time delays and probe pulse 
widths and integrating the tr- ARPES signal over energy. 

We illustrate the use of the TDMD by demonstrat- 
ing how it can be used to map out the unoccupied dis- 
persion in a pump-probe setup on a model system with 
electron-lattice interactions appropriate to a large class 
of complex materials. The TDMD exhibits two ubiqui- 
tous phenomena in time-resolved experiments; both are 
due to electron-phonon interactions. First, the phonons 
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FIG. 1. a) (Color online) Equilibrium TDMD n^(t) far before the pump arrives. The gray area indicates where quasiparticle 
states are occupied. The Fermi surface is sharply defined for the parameters presented (g 2 — V0.021/ nn , T — 0.04y nn ,Q = 
OAVnn)- b) Change in the momentum distribution from equilibrium Sn^t) after a pump at t ~ 4.8r p . The arrow indicates the 
direction of the applied electric field. 



dissipate the energy delivered to the electronic system by 
the pump. Second, after pumping, oscillations related to 
the phonon frequencies (and independent of the initial 
pump parameters) are commonly observed in the result- 
ing spectra measured by the probe [7HT0]. We will show 
that both of these features are clearly seen in the TDMD. 

We study a two-dimensional system of electrons cou- 
pled to a bath of non-dispersive phonons with frequency 
Q via a constant coupling g (known as the Holstein 
model), HQ. 
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where c k (b q ) annihilates an electron (phonon) of mo- 
mentum k (q). 

While the TDMD is useful independent of the under- 
lying model, here we use a band structure motivated 
by the transition-metal oxides: a 2D tight-binding e(k) 
with (next) nearest-neighbor hoppings V nn and V nnn = 
0.3V nn , and \i = — 1.02V^ n . We use the conventions 
h = e = c = 1. Below, we will report timescales in units 
of a characteristic phonon timescale r p = I/O. For exam- 
ple, a phonon with an energy of O = 10 meV has r p ~ 0.4 
ps. Similarly, there is a characteristic electron timescale 
r e = 1/V nn ~ 16 fs for V nn = 250 meV. Although in this 
work we will explicitly examine higher energy phonons 
for numerical stability reasons, the results below occur 



due to the relevant energy scales in the problem, and 
will simply be rescaled to lower phonon energies. 

To describe the non-equilibrium pump-probe process, 
we propagate the system on the Kadanoff-Baym-Keldysh 
(Keldysh) contour [12], which has been described in detail 
elsewhere [T3HT5] . The non-equilibrium Keldysh formal- 
ism has been successfully used to describe the time evolu- 
tion of correlated electrons within dynamical mean-field 
theory [16-19 . The electric field is included through the 
standard Peierls' substitution k — » k — A(t), where the 
vector potential A is related to the applied electric field 
F via —dA(t)/dt — F and we work in the Hamiltonian 
gauge. 

Within the Migdal limit, we consider the renormaliza- 
tion of the electron Green's function by a single phonon 
emission/absorption and ignore the phonon self-energy 
(as is done in conventional Migdal-Eliashberg theory). 
In this limit, we solve the Dyson equation 



G k (t, t f ) = Gl(t, t f ) + j dhdt 2 Gl(t, OEfa, t 2 )G, 
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with the self-energy E(t,t') = W 2 Ek D °(t> O^kfo 0> 
where G^(t,t f ) and D°(t,t') are the non-interacting elec- 
tron and phonon Green's functions, respectively [T3j . 
The Dyson equation can be solved by recasting it as 
a matrix equation [15 or by decomposing into Volterra- 
type equations through the Langreth rules [20 , which is 
the approach used here. 

Figure [l] shows the gauge- invariant TDMD n^t) = 
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iG k + A( t) (M ' ^ VMfor g = VomV nn , O = OAV nn 
and at a temperature T = 0.04V nn . Fig. [l^i) shows 
nk(t) at times far before the pump, where the system 
is fully described by the equilibrium problem. To clearly 
show the changes for the small pump fluences considered 
(Fmax < 4Kn) Fig. [!]□) shows the change in the TDMD 
from equilibrium: Sn^t) = n^(t) — n^(t — »> — oo). 

We apply a pulse centered at time t = in 
the (11) (diagonal) direction of the form A(t) = 
(^max/^) exp(— t 2 / (2a 2 )) sin(cj^), with a maximum 
field strength of F max = AV nni a frequency cj^ = 2T/ nn 
and <j w 0.4r p . At the time shown in Fig.[l]b), the pump 
has passed, and the system is relaxing towards equilib- 
rium. The C4 symmetry of n^(t) in equilibrium is broken 
by the field (which points in the (11) direction as shown 
by the black arrow), and this is reflected in the transient 
change. In fact, this can be used to emphasize particular 
regions of interest in the Brillouin zone. By rotating the 
field, the regions exhibiting the largest change in n^(t) 
will shift with the field direction, providing access to new 
regions by emphasizing the experimental signal near dif- 
ferent momenta. 
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FIG. 2. a) (Color online) Change in the TDMD n k (t) from 
the equilibrium value (at t — 0) for the points in the Brillouin 
zone indicated in b) . Symbols are calculated data points (for 
clarity, not all are shown) , solid lines are fits to an exponential 
decay using a fixed relaxation rate 1/rk = — 21m S(cj = ek). 
c) Real and imaginary parts of the equilibrium self-energy 
S(cj) for the parameters in a) (g = V0.02V^ n , T = 0.04V nn , 

n = 0AV nn ). 



First, we consider the anisotropic redistribution of 
quasiparticles as a function of time. In Figure^), we 
show 5n\z(t) at the momenta indicated in the inset; the 
blue (black) momenta have the same energy ek but due 
to the C4 symmetry breaking have different time traces. 
From the figure, one can clearly see that there is a dif- 
ference in the overall magnitude of the change at the 



equivalent momenta (with the same color). The solid 
lines are fits to the amplitude of a simple decaying ex- 
ponential, where we have taken the decay constant at a 
momentum k from the on-shell equilibrium self-energy 
l/rk(cj) = — 21m E(co> = ek) shown in Fig. J2Jd), as sug- 
gested by a recent study [2 1 . The remarkable agreement 
with the calculated data points confirms that it is in- 
deed the equilibrium self-energy that controls the decay 
rate. Since the equilibrium self-energy is C4 symmetric, 
the anisotropic redistribution in Sn^t) is caused by the 
C4 symmetry breaking of the applied pump field. 
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FIG. 3. a) (Color online) Oscillatory part of 5rik(t) for k 
along (11) at the energies ek indicated (offset for clarity), b) 
Fourier transform of a) . Vertical lines indicate the frequencies 
corresponding to the strong oscillations on the unoccupied 
side, u — ek + H and uj — ek — W- + In both a) and b) the 
oscillations on the unoccupied side close to the Fermi level 
(blue curves) have been scaled down by 4 for visibility, (see 
text for details). 



A ubiquitous feature observed in pump-probe experi- 
ments is oscillations [7 -10 of the time traces dominated 
by a few characteristic frequencies. These are clearly ob- 
served as rings in Sn^t) shown in Fig.JT] as well as oscil- 
lations in the time traces in Fig. [2] We further emphasize 
these in Fig. |3j where we have subtracted the decaying 
exponential from several points along the Brillouin zone 
diagonal in Fig. [l] From Fig. [3^t) , one can observe two 
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main features. First, the oscillation frequencies in the 
time traces depend strongly on momentum. Second, al- 
though there is only a single phonon frequency in our 
model, several oscillation frequencies can be observed. 
Figure [3]o) shows the Fourier transform power spectrum 
of the oscillations in Fig. [3^). Each curve has readily vis- 
ible maxima, in addition to some smaller structure across 
the frequency spectrum. For each energy (momentum) 
ek, the vertical lines indicate the energies uj = €k + Q and 
uj = ek — W- + ^, where W is the energy at the bottom 
of the band (at k = V). 

The temporal dynamics of the Holstein model leads to 
this rich physics. The decay and the oscillations observed 
are momentum (energy) dependent, and their amplitudes 
depend on the electron-phonon coupling, as well as the 
applied field. Describing the full temporal dynamics of 
n^(t) is complicated, yet we can gain some insight by 
focusing on the times where the pump is off, while the 
system is still displaced from its equilibrium configura- 
tion. The subsequent relaxation occurs from the new 
configuration through the equations of motion. In equi- 
librium, the TDMD is balanced by equal rates of in- and 
out-scattering; once the pump breaks this symmetry, the 
system relaxes back to the equilibrium according to the 
dynamics contained within the self-energy and Green's 
functions. In addition to the relaxation, the TDMD os- 
cillates at the band energy shifted by the characteristic 
energies in the problem (the phonon frequency and band- 
width) at which the self-energy is largest. This is shown 
analytically at T = in the supplemental information; 
the full calculation requires the numerical methods used 
here. 

At low temperatures, the real part of the self-energy 
Re E is large and peaked near four frequencies, uj — ±fi 
and uj = W± ± f2, where W± is the upper (lower) band 
edge (see Fig. |2j. The resulting n^(t) oscillates at the 
band energy ek shifted by those frequencies. In princi- 
ple, all four frequencies should be observable in the power 
spectrum of Fig. [3j however, two frequencies (with the 
positive sign above) correspond to phonon absorption, 
which is small at low temperatures, and thus the corre- 
sponding peaks in the power spectrum are reduced sig- 
nificantly. The two strong peaks expected for the un- 
occupied side are indicated in Fig. [3] by vertical lines, 
which agree well with the observed frequencies, up to a 
small shift in energy due to Re E. An interesting con- 
sequence of the expected oscillation frequencies is that 
if one considers a measurement near the Fermi level, the 
oscillations will appear to be strongest at just the phonon 
frequency (as observed in tr- ARPES [7HT0] ) . 

The simple dependence of the oscillation frequencies on 
the underlying band structure suggests a novel method 
to measure the electron dispersion by looking at the oscil- 
lations in the TDMD. The oscillations are most clearly 
visible in the unoccupied region of the Brillouin zone, 
which is exactly the region that traditional methods for 



measuring the dispersion have difficulty accessing. By 
tracking the oscillation frequencies as a function of mo- 
mentum one can directly map out the dispersion. In 
Fig. [4j we plot the Fourier transform power spectrum in 
the unoccupied portion of the Brillouin zone along the 
zone diagonal and zone face. This shows most clearly 
the two strong peaks in the Fourier transform associ- 
ated with phonon emission, although a weaker line cor- 
responding to phonon absorption is also visible below the 
two main lines. All three lines, up to the constant shifts 
from the real part of the self-energy described above, fol- 
low the unoccupied dispersion. Thus, from the maxima 
in the power spectrum the unoccupied dispersion can be 
directly measured. 

The above analysis, applied to theoretical results, can 
be applied to experimental results equally well, although 
one must take care to include the a 2 F(Q) for the ma- 
terial under study. Time-resolved Compton scattering 
also requires the introduction of the scattering matrix 
elements, similar to what has been addressed in equi- 
librium Compton scattering [3 . Additionally, Compton 
scattering measures a single momentum direction, aver- 
aging over the directions not selected by the field. This 
allows for selectivity in the direction of interest, focus- 
ing on the regions of the Brillouin zone that one wishes 
to measure. Finally, there is a distinction between the 
crystal momentum in the TDMD and the real momen- 
tum measured by Compton scattering. The inversion 
from real momentum to crystal momentum has been ad- 
dressed in some detail previously. [22] 

Experimentally, the width of the probe pulse limits 
the temporal resolution of the TDMD. Thus the Fourier 
transform to frequencies is limited in turn to several hun- 
dred meV around the Fermi level. However, in many sys- 
tems of interest, the features one wishes to investigate lie 
near the Fermi level, both above and below. For example, 
in the high-Tc cuprates, several features are predicted to 
lie within this energy range. [23] In other gapped systems, 
without the knowledge of the spectra above the Fermi 
level it is impossible to correctly assign the gap magni- 
tude. Similar effects due to spin fluctuations, phonons, 
and other interesting collective and bosonic features lie 
at low energy, and are thus accessible by this technique. 
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FIG. 4. (Color online) Fourier power spectrum of the oscillation frequencies along the T-M and M-X directions in the unoccupied 
Brillouin zone. The two strong peaks in the Fourier spectra follow the lines corresponding to phonon emission processes; a 
small tertiary line of peaks can be seen at lower frequencies corresponding to a phonon absorption process. The Fermi surface 
is shown in green on the base plane. 
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